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We analyse the redshift-space (z-space) distortions of QSO clustering in the 2dF 
QSO Redshift Survey (2QZ). To interpret the z-space correlation function, £(er, 7r), we 
require an accurate model for the QSO real-space correlation function, £ (r) . Although 
, a single power-law £(r) oc r~ 7 model fits the projected correlation function (w p (a)) 

at small scales, it implies somewhat too shallow a slope for both w p (a) and the z- 
space correlation function, £(s), at larger scales (> 20 ft, Mpc). Motivated by the 
form for £(r) seen in the 2dF Galaxy Redshift Survey (2dFGRS) and in standard 
ACDM predictions, we use a double power-law model for £(r) which gives a good fit 
\ to £(s) and w p (a). The model is parametrized by a slope of 7 = 1.45 for 1 < r < 

10/i _1 Mpc and 7 = 2.30 for 10 < r < 40/i~ 1 Mpc. As found for 2dFGRS, the value 
of /3 determined from the ratio of £(s)/£(r) depends sensitively on the form of £(r) 
assumed. With our double power-law form for £(r), we measure P(z = 1.4) = 0.32±g;??. 
Assuming the same model for £(r) we then analyse the z-space distortions in the 2QZ 
£(cr, 7r) and put constraints on the values of and (3{z = 1.4), using an improved 



Oh 



version of the method of Hoyle et al. The constraints we derive are fi^i = 0.35^q 'jf, 
(3(z = 1.4) = 0.50tai5, in agreement with our t;(s)/t;(r) results at the ~ 1 a level. 
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1 INTRODUCTION 

Understanding the formation, clustering and evolution of 
large-scale structure (LSS) has been one of the major quests 
in the history of cosmology. One of the ways of quan- 
tifying the inhomogeneity of the Universe on the largest 
scales is to count the number of pairs of galaxies or quasi- 
stellar objects (QSOs), as a function of separation between 
them. The resulting correlation functions from such statisti- 
cal analyses a ppear to show a power-law behaviour for scales 
<10/i _1 Mpc JPeebleslll974lh 

Statistical tools such as correlation functions or power- 
spectra have been used in many cosmological studies. Un- 
fortunately, the underlying dark matter clustering may dif- 



fer from that of the luminous matter, a difference that can 
be quantified by the bias (Kaiser 1984). Black-hole models 
make a wide variety of predictions for the amplitude and 
form o f the bias (e.g. Kau ffmann & Haehnelt 2000). Else- 
where, ICroom et all (I2005T) have tried to constrain models 
for the QSO bias via studies of QSO clustering evolution. In 
this paper, we shall use linear dynamical models to inves- 
tigate QSO clustering in "redshift-space", with the aim of 
constraining the amplitude of the QSO linear bias (6) and 
the value of the cosmological density parameter fi^. 

To compute the distance to the QSOs one needs to as- 
sume that their observed redshifts are only due to the expan- 
sion of the Universe. However, the observed redshifts also in- 
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elude contributions from the peculiar motions of the QSOs. 
Therefore, the estimated distances are said to be measured 
in redshift- space, rather than in real-space. Peculiar veloci- 
ties introduce distortions in the measured clustering pattern. 
If a given ensemble of QSOs has, on average, a spherically 
symmetric clustering pattern in real space, but the QSOs 
have a very large velocity dispersion, then the clustering 
signal measured in z-space will be smeared along the line- 
of-sight. These features are commonly referred to as "fingers- 
of-God". At large scales, gravitational instabilities can also 
generate peculiar motions due to a coherent infall of QSOs 
into the potential well of higher density regions, causing a 
flattening of the clustering in the z-direction. The presence 
of the these dynamical distortions, namely the large-scale 
flattening at large scales, has been observed in galaxy sur- 
veys such as the 2dF Galaxy Redshift Survey (2dFGRS, 
Hawkins et al, 2003). 

Besides the dynamical distortions, geometric distortions 
in the .z-space clustering pattern also occur if a wrong cos- 
mology is assumed to convert the measured redshifts into 
comoving distances. This is because the separations esti- 
mated along and across the line-of-sight have different de- 
pendences on cosmology, and hence they will change dif- 
ferently if incorrect cosmologies are assumed. The use of 
geometric distortion s to st udy cosmology was pioneered by 
lAlcock fe Paczvnskil dl979l) . who demonstrated that they are 
potentiall y a powerful cosmological test for a non-zero A. 
However, iBallinger. Peacock fc Heavens! il996T) noted that 
the geometric and dynamical effects may be degenerate. 

To break this degenerate constraint between A and bias, 
complementary constraints on these pa rameters can be con - 
si dered simultan e ously. As shown by iHovle et alJ J2002fl . 
or lOutram et al] (I2004T) . orthogonal constraints can be ob- 
tained from linear evolution theory of cosmological density 
perturbations. In this paper we develop a similar analysis 
using the full catalogue of the 2dF QSO Redshift Survey 
(hereafter 2QZ). We aim to obtain constraints on the QSO 
bias at z « 1.4 and on the value of the matter density of the 
Universe from a detailed analysis of the z-space distortions 
in the 2QZ clustering, combined with bias constraints from 
linear theory. 

We find that, to obtain an adequate model describing 
the clustering along and across the line of sight, we need an 
accurate description of the real-space correlation function, 
£(r). We develop a model for £(r) based on the results from 
the 2QZ z-space and projected correlation functions and use 
it as an input to the two-dimensional z:-space distortions 
model. 

The structure of this paper is as follows: in Section 2 
we introduce the 2QZ sample and outline some of the most 
relevant aspects of the survey. In Sections 3 and 4 we start 
by introducing our correlation function measurements, and 
derive £(r) models suitable for explaining our observations. 
We then compare our results with the clustering analysis of 
the 2dFGRS survey (Section 5) and with predictions from 
CDM non- linear clustering models (Section 6). In Section 
7 we discuss further constraints that can be obtained on 
/3(z = 1.4) and from fitting the z-space distortions in 
the two-dimensional z-space corre lation function. W e also 
compare our method with that of IHovle' et al] (120021) . m a 
similar study. Finally, in Section 8, we draw the conclusions 
of this work. 



Throughout this paper we use two different cosmologies 
to compute the comoving distance as a function of redshift: 
an Einstein-de Sitter Q,^ = 1.0 cosmology (EdS), and a flat, 



0.3, Q\ = 0.7 cosmology (A). 



2 THE DATA 

The 2QZ has presented the scientific community with a 
multitude of new and interesting results. Comprising 23338 
QSOs with a limiting magnitude of bj = 20.85, it is the 
largest homogeneous completed QSO survey by more than 
a factor of 50, at such faint magn itudes and with such high 
surface density (35 QSOs deg" 2 ) l|Croom et al]l2004aT ). The 
high number of QSOs and their spatial density make the 
2QZ a powerful tool for studying QSO clustering, in order 
to determine the real space clustering and model the z-space 
distortions in the clustering pattern. 

The survey area consists of two declination strips, 75° x 
5° each. One is centred on 5 = —30° and extends from a = 
21 h 40 m to q = 3 h 15 m and the other is centred on S = 0° 
and extends from a = 9 h 50 m to a = 14 h 50 m . (Note that 
these are values in the B1950 coordinate system.) 

The QSOs have 18.25 < &.j < 20.85 and were se- 
lected using the ultra-violet excess technique in the plane 
u — bj : bj — r. Due to the colour selection used, they lie 
in the redshift range 0.3< z< 3.0. The number density of se- 
lected QSOs decreases abruptly close to both redshift limits. 
At low redshifts, two effects contribute to this decrease: 1) 
the host galaxies' colours contaminate significantly the ob- 
served QSO colours, since faint QSOs are being detected; 2) 
the host galaxy appears extended and these objects will be 
missed, since only point-sources are considered in the sur- 
vey. The high redshift limit is due to the Lyman forest that, 
at those redshifts, lies in the optical and therefore causes 
a significant absorption in the u band. The analysis per- 
formed here only includes QSOs in the range 0.3 < z < 2.2, 
where the photometr ic completeness of the survey is > 90% 
l|Croom et aljl2004al) . 

Around 47% of the colour - selected candidates were 
spectroscopically identified as QSOs, with the criterion 
being the presence of broad (> 1000 kms^ 1 ) emission 
lines. The remaining objects are either stars, white dwarfs, 
narro w-line galaxies, cataclysmic variables ( Marsh et al] 
2002), or their identification was not possible from the spec- 
tra. 

The spectra were obtained using the 2-degree field in- 
strument (2dF) at the Anglo- Australian Telescope (A AT). 
The design of this multi-fibre spectrograph allows the ac- 
quisition of 400 spec tra simultaneousl y over a 2° diameter 
circular field of view l|Lewis et al.l20 02V To cover the survey 
area with the 2dF field, a tiling algorithm that maximizes 
the number of galaxies and QSOs per pointing was used. The 
tiling of the 2dF fields has a significant effect on the angular 
selection function of the survey. A very detailed analysis of 
the angular and radi al QSO selection fun ction s used in this 
work ca n be found inJOutran^^L J2003J) and lCroom et all 
l|2004al) . Following Iputram et al] J2003I) . we introduce the 
redshift cuts 0.3 < z < 2.2, which allow us to consider clus- 
tering statistics without the need for a weighting scheme. As 
a consequence, the final QSO catalogue used in this study 
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comprises 19549 QSOs (8704 and 10845 QSOs in each of the 
strips) . 



3 THE REDSHIFT-SPACE TWO-POINT 
CORRELATION FUNCTION 

3.1 Estimation of The Correlation Function 

Statistics such as the two-point correlation-function (£) can 
be used to describe the clustering and the distribution of 
QSOs. f(a;) is the excess probability (due to clustering) of 
finding a pair of objects at a given sepa ration - x - re lative 
to a Poisson (unclustered) distribution (Peebles 1980). 



SP = n8V{l + £(x)) 



(1) 



is the probability that one has of finding a QSO, at separa- 
tion x from a neighbour, located in a volume SV. n is the 
mean number density of QSOs. 

To evaluate how much the QSO distribution deviates 
from an unclustered set, a random ensemble of points prob- 
ing the same volume as the QSOs is generated. This random 
distribution must describe all the characteristics of the QSO 
ensemble, except for its clustering, and therefore it must be 
generated taking into account the sa me selection function s 
as the QSO catalogue. For details see lCroom et al.l (|2p04a). 

To evaluate £ different estima tors have been de veloped. 
Here we use the one developed bv lHamiltonl (ll993T> : 



< DD(x) >< RR(x) > 
< DR(x) > 2 



(2) 



where < DD(x) > is the number of QSO-QSO (double 
counted) pairs separated by x, < RR(x) > the number of 
random-random pairs separated by x and < DR(x) > the 
number of QSO-random pairs also separated by x. 

At small and intermediate separations, where the pairs 
in each bin are independent, the errors on the correlation 
function are determined using the Poisson estimate. Their 
computation is quite straightforward: 



A£ = (1 + 



<DD> 



(3) 



The factor \/2 arises from the fact that each pair is counted 
twice. At large scales, the pairs in each separation bin 
are not independent, and hence Poisson errors are not 
applicable and undere stimate the true errors. Following 
Shan ks~fc Bovld il994h . at large scales the errors can be 
computed by A£ = (1 + vV^Qi where Nq is the total 
number of QSOs used for computing £. We use this estimate 
at separation bins where the total number of pairs exceeds 
the number of QSOs in the survey. Alternatively, the errors 
could be estimated by considering the dispersion in either 
mock QSO catalogues derived from iV-body simulations, or 
subsamples of our QSO ensemble. However, as shown by 
iHovld (2000), the Poisson and Nq errors, as used here, are 
in good agreement with the errors obtained with these other 
estimators. Note that the correlation function estimates may 
be correlated bin-to-bin. 
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Figure 1. The z-space two point correlation function £(s) for 
the complete 2QZ survey, measured with the Hamilton estimator, 
assuming A. The solid line represents the best fitting power-law 
for 1 < s < 100?t _1 Mpc. The dashed line is the best fitting 
power-law for 1 < s < 25/i _1 Mpc. 



3.2 The 1-D z-Space Correlation Function £(s) 

A detail ed description and a nalysis of the 2QZ £(s) can be 
found in lCroom et~ail (|2005ll . They find that the 2QZ £(s) 
is not well described by a single power-law model. Fig. 
shows the 2QZ £(s) measured using the Hamilton estimator, 
assuming the A cosmology to compute the r — z relation. 
The lines represent the bes t fitting power-law models for 
different scales, computed bv lCroom et alJ <l2005l) . The best- 
fitting power-law over scales 1 < s < 100ft _1 Mpc 1 , £(s) = 
(s/5.55) -1 ' 63 , is shown by the solid line. The dashed line 
represents the best fit at scales 1 < s < 25/i _1 Mpc; £(s) = 
(s/5.48)" 1 - 20 . 

A noticeable feature in Fig. Q quantified by the dis- 
crepancy between the two best fitting power-laws, is the 
non-power law shape of £(s), being flatter at small scales 
relative to the slope at larger separations. This small-scale 
flattening is frequently observed in redshift surveys and is 
usually attributed to the effect of z-space distortions, namely 
as a result of the QSOs' velo city dispersion in the ir local rest- 
frame. As demonstrated bv lCroom et alJ (l2005f) . this small 
scale flattening is not due to the fibre collisions of the 2dF 
instrument. In addition, it was also found that the major- 
ity of QSO pairs at small angular separations actually have 
quite large radial separations, hence not affecting the shape 
of £(s) at small scales. We found that the correlation be- 
tween the £(s) bins was small. The covariance matrix shows 



1 In this paper, we will denote all the separations measured in 
redshift space with s and the ones measured in real space with r. 
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that the correlation between adjacent bins is measured to 
be only ~ 10% over the range of scales fitted (e.g. Myers et 
al. 2005). 



4 A COHERENT PICTURE OF THE QSO 
CLUSTERING 

Studying the real-space clustering of the QSOs is, in itself, 
a subject of obvious interest: it allows the direct study of 
QSO clustering and its bias relative to the underlying dark 
matter distribution, providing a picture of how the forma- 
tion and evolution of QSOs take place. Nevertheless, our 
motivation for studying f(r) also includes the need for an 
adequate and trustworthy amplitude input for modeling the 
shape of the clustering in orthogonal directions. Comparing 
the distortions in the clustering measured along and across 
the line-of-sight will allow the measurement of cosmological 
and z-space distortion parameters. 

Let the separations along and across the line of sight be 
defined as: 



s 2 



si 



(Sl + s 2 ) 



(4) 
(5) 



where si and s 2 are the distances to two different QSOs, 
measured in z-space, and 9 the angular separation between 
them. 

Since the effects of z-space distortions are purely radial, 
the real-space clustering can be inferred from the projection 
of £(<7, 7r) along the a direction 2 . The projected correlation 
functio n is obtained b y integrating £(<r, 7r) along the n di- 
rection JPeebleslll98rjl) : 

w p {o) =2 Z(a,n)dn (6) 
J o 

When implementing this integration it must be trun- 
cated at some value n cu t, where £ becomes negligible. If very 
large scales are included, the signal amplitude will become 
dominated by noise. On the other hand, if the integral is only 
performed at very small scales, then the projected correla- 
tion function will be systematically underestimated. It was 
found that, in the present survey, the results did not seem 
sensitive to ir cu t values greater than 70 /i _1 Mpc, which was 
the value used. Therefore: 



w p (a) = 2 



i;(o", n)dir 



(7) 



Fig. |21 shows w p (a)/a measured from the 2QZ cata- 
logue. The error bars represent Poisson errors. 

Since w p (a) describes the clustering in real-space, the 
integral in equation |S| can be written in terms of the real- 
space correlation function £(r) foavis fc Peebleslll985t) : 



w p (a) 



\f>' 



zdr 



(8) 



Substituting £(r) for a power law of the form [j-j 
equation |H1 w p (a) will be given by: 



£(<7, n) will be considered in detail in Section 7. 







n° m = i-o 


Q° m = 0.3 


T'Q (h 


_1 Mpc) 


o w7~|"0.33 
°- 1 1 -OAT 


4 Q6+ - 54 
^• MD -0.56 


7 




2 09+ ' 21 


i or-+0.13 

1,8& -o.io 


X 2 ■ 

V HI 1 It 


(reduced) 


1.61 


0.93 


d.o.f. 




6 


6 



Table 1. The values of ro and 7 from the best fitting power-law 
model to w p (cr)/cr, for 6 degrees of freedom (d.o.f.), corresponding 
to the solid lines on Fig.|5] If the EdS model is assumed to convert 
the redshifts to distances, the fit is performed between 1 h _1 Mpc 
and 25 fc~ 1 Mpc. Assuming the A cosmology the range 1.2 < a < 
30/i Mpc was taken. 



(*) = 



'rQ)r(^)' 



r i 



(9) 



where F(x) is the Gamma function computed at x. 

Therefore, if the real-space correlation function £(r) is 
well approximated by a power-law, then its slope will be the 
same as in w p (a)/a, and its amplitude proportional to the 
w p (a)/a amplitude. Fig.|!2]also shows the best fitting power- 
law to w p (a)/a (solid line). The fitting range is scaled ac- 
cording to the cosmology assumed. The range used is 1 < 
a < 25/i _1 Mpc, assuming EdS, and 1.2 < a < 30h -1 Mpc, 
assuming A. The corresponding best fitting values of ro and 
7 are summarized in Tabled The errors on these parame- 
ters can be inferred from the confidence levels in the [ro,7] 
plane, shown in Fig. [3] The errors in Tableland Fig^may 
be underestimates because of the possibility that the w p (a) 
points are correlated. However, in jacknife tests the corre- 
lation between adjacent bins is measured to be only ~ 20% 
over the range of scales fitt ed. The corr elation analysis of 
Hubble Volume simulations l)Hovlell2000h also suggest that 
the w p (a) points are approximately independent. 

The close-to-unit values of the reduced x indicate that 
the power-law prediction for £(r) corresponds to a rea- 
sonable description of w p (a) within the limited range it 
was fitted. If EdS is assumed, then the best fitting £(r) 
power-law model has the form: (r/3.77) -2 ' 09 , while if A 
is assumed instead, the best fitting power-law model is: 
£(r) = (r/4.96)" 1 - 85 . There are deviations from the best- 
fitting power-law model at scales larger than the ones used 
in the fit; the w p (a) clustering amplitude is lower than pre- 
dicted from the best-fitting model, for er> 30 ft -1 Mpc. 



4.1 Does a single power-law £(r) model explain 
both the w p (a) and £(s) results? 

The simplest explanation to account for the different slope 
of £(s) at 1< s< 10 h~ Mpc, compared to the slope at larger 
scales, is the effect of z-space distortions. The small-scale 
random motions of the QSOs lead to a deficit of clustering 
amplitude measured at scales <10/i _1 Mpc. In addition to 
this effect, the £(s) clustering amplitude on large scales is 
affected by the infall of QSOs into the overdense regions. 

According to linear th eory, the rela tion between £(s) 
and £(r) is simply given by feaiseilll987ft : 
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Figure 2. The projected correlation function measured with the 
Hamilton estimator and the best fitting power-law. The top panel 
shows the case where EdS was assumed, and the bottom panel 
the result assuming A. The solid lines represent the best fitting 
power-laws, that correspond to £(r) = (r/3.77) — 2 09 assuming 
EdS and £(r) = (r/4.96) -185 assuming A. 




(10) 



where f3 w 0.™ /b. 

Therefore, the ratio between the z- and real-space cor- 
relation functions gives an estimate of the infall parameter 
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Figure 3. The best fitting values of tq and 7 and the respective 
two-parameter confidence levels obtained when fitting a model of 

to w p (<j). The solid contour 



the form r^a 1 7 



lines represent the results from fitting the function obtained as- 
suming A, and the dotted line for the EdS cosmology. The best 
fitting pairs of values [ro ,7] are indicated with X . 



j3(z = 1.4) which, for a given cosmology, allows the deter- 
mination of the bias b. 

With the current data, it is not trivial to find the real- 
space correlation function without making some assump- 
tio ns. Inversion method s o f w v (a) like the ones described 
bv lSaunders et alJ <ll992l) or lRatcliffel {596) are not feasible 
with this survey, since the errors associated with f (r) would 
be extremely large. 

If we approximate £(r) by a single power-law, then the 
best fitting power law w p (a) model can be used to constrain 
/3(z — 1.4). This can be done by computing the ratio be- 
tween the measured values for £(s) and the values our £(r) 
model takes, at large scales. Fig.0|shows the ratio £(s)/£(r) 
as a function of separation. For a given value of the sepa- 
ration, the error on £(s) can be quantified using the Pois- 
son estimate. The uncertainty in the £(r) model is deter- 
mined from that in the parameters ro and 7. However, as 
these parameters are correlated (see contours in Fig.[3Jl, the 
value found underestimates the true uncertainty. We there- 
fore rescale the error in f(r) by comparison with the £(s) 
Poisson errors. To estimate the error in the ratio £( s )/£( r )i 
we then add in quadrature the error in £(r) and the Poisson 
errors in £(s). The dashed line in Fig. 2] shows the best fit- 
ting £(s)/£(r) value at scales larger then 10 /i _1 Mpc, and 
the shaded region is the 1 a confidence interval in the fit. 
The best fit produces (3(z = 1.4) = 0.871^?. At scales 
>10/i _1 Mpc, non-linear effects due to peculiar velocities 
have a negligible effect and the z-space contribution comes 
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Figure 4. £(s)/£(r) measured from the 2QZ survey, taking the 
measured f (s) values and the single power-law £(r) model derived, 
assuming A. See text for a full description of how the errors are 
computed. The fit to the function was performed on scales where 
the 2-space distortions are only affected by the large-scale infall 
and are not contaminated by random peculiar motions of the 
QSOs. The dashed line and the shaded region represent this best 
fit and its 1 a confidence level. This fit corresponds to a value of 
P(z = 1.4) = 0.87 " 
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only from the gravitational infall quantified by /3(z). This 
value of /3(z = 1.4) is larger than other r esult s also from the 
2QZ s urvey, found bv lHovle et al.l i2002T) and lOutram et ail 
(2004). We should point out that, using single power-law 
fits to the w p (a) result of the 2dF Galaxy Redshift Survey 
(2dFGRS; Hawkins et al. 2003) also gives high values of (3 
from computing £( s )/£( r )- 

A model for £(s) can then be constructed by decom- 
posing the separation in a and 7r and adding the z-space 
distortions in tv). Then, in order to obtain £(s), £(cr, 7r) 
is averaged in annuli of constant s = \f a 2 + ty 2 . The model 
used to compute £(cr, n) is discussed in Section 17.31 Fig. Q>1 
shows the comparison of £(s) models with the 2QZ £(s). The 
circles in the plot show the 2QZ £(s), assuming the A cos- 
mology. The dotted line shows the input £(r) model, which is 
the best-fitting power-law to the w p (a). Note that this model 
lies well below the observed £(s) at all values s> 10 /i _1 Mpc. 
The dashed line represents that same power-law, scaled us- 
ing equation 1101 taking the value of (3(z) — 0.87 to quantify 
the effect of the linear z-space distortions. Including the dis- 
tortions due to small scale pairwise velocity dispersion of the 
QSOs, will lead to the solid or the dash-dotted lines. The 
solid line shows the effect of small-scale peculiar velocity dis- 
persion < vol y 1 ^ 2 — 800 kms" 1 . This value is dominated 
by the rms pairwise redshift-error of 600 kms" 1 measured 
from repeated observations with the remainder allowing for 
the intrinsic velocity dispersion of the QSOs. The effect of 
a larger velocity dispersion of 1500 kms -1 is shown as the 
dot-dashed line. 

It can be seen that the lower velocity dispersion of 
800 kms -1 gives an adequate fit to the data in the 
s< 10 /i _1 Mpc region; the tendency of this model to over- 
estimate the points at 2<s<4 ft -1 Mpc is not resolved by 
moving to a bigger velocity dispersion which is seen to de- 
grade the fit at s ~ 10 ft _1 Mpc. Although the \ 2 ls ac ~ 
ceptable, the slope of £(s) at s> 10 /i _1 Mpc is steeper 
than the single power- law model derived from w p (a). It 
proved impossible to increase the slope in £(r) in this range 



Figure 5. The circles show the measured £(s) from the 2QZ, 
assuming A. The dotted line is the best fit single power-law £(r) 
model to w p (o). The dashed line is the same model scaled to 
account for the linear bias {/3(z) = 0.87). The solid and dash- 
dotted lines also include the small-scale peculiar velocities with 
velocity dispersion of the order < w 2 > 1 / 2 = 800 kms -1 and < 
w 2 > 1 / 2 = 1500 kms -1 , respectively. 



while maintaining an acceptable fit to £(s). This tendency 
of single-power law £(r) models to give poor fits to £(s) 
has been seen i n previous galaxy red shift surveys such as 
Durham/UKST Ratcliffe^jtaljjm)98|) and more recently in 
the 2dFGRS lHaw kins et al] fcOOSTT where a non-power-law 
shoulder is seen rising above the small-scale power-law at 
s ~ 10 /i _1 Mpc, before steepening again at larger scales. 



4.2 A double power-law £(r) model: another 
explanation of the w p (a) and £(s) results 

lHawkins et al] i2005t) presented the correlation function 
analyses of the 2dFGRS, both in real- and z-space. Including 
data from 166000 galaxies with a mean redshift of z m 0.11, 
this study has a statistical weight significantly higher than 
the 2QZ. By inverting the projected correlation function, 
they found that the slope of the real space correlation func- 
tion of the 2dF galaxies varies with scale. In order to quan- 
tify the change in slope as a function of scale, we re-fitted 
their £(r) data. At scales < 1 /i _1 Mpc the shape is si milar to 
the c ommonly observed 7 w 1.8 small-scales slope dPeeblesI 
Il974l) however, at intermediate scales (1 — 10 /i _1 Mpc), that 
are accessible with the 2QZ survey, the £(r) data is fitted 
by a shallower 7 = 1.45 power- law, while at larger scales the 
correlation function is again steeper. Motivated by these re- 
sults we now assess if the addition of a break in the £(r) 
shape at ~ 10 h~ Mpc can still provide a good description 
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of the QSO w p (a) and { (s) results, whilst accounting for the 
issues discussed in the previous section. 

We thus perform a new parametric fit to w p (a), consid- 
ering a double power-law model with a break at 10/t _1 Mpc 
(for the A cosmology). Then, for a grid of parameters that 
characterize the slope and amplitude of the power-laws, we 
project the functions describing { (r) using equation|H|to get 
the respective w p (a). The best fitting model is found by 
performing a \ 2 fit- The results are shown in Fig. |S| The 
value used for the upper limit of the integral in equation [S] 
was computed from the n cut value used in equation |7| , to 
determine w p (a) from the data. 

The top plot shows the best fitting £(r) model to w p (a) 
(solid line). The parameters that describe the two power 
laws are: for scales r < 10 /i _1 Mpc, r = 6-0t°.l /i _1 Mpc 
and 7 = 1A5±° 2 2 7 T ; for r > 10 h^Upc, r = 7.25 /i _1 Mpc 
and 7 = 2.30^q'J3. The fit is performed in the range 1.0 < 
a < 40/i _1 Mpc. We fit the slope and amplitude of the two 
power-laws, keeping the break between the two-power laws 
at 10/i _1 Mpc and ensuring that the function is continuous 
across the break. The reduced Xmin °f this fit is 0.89 (6 
d.o.f.). The bottom plot is the measured 2QZ w p (a)/a over- 
plotted on the projection of the best fitting f (r) model (solid 
line) . 

This model is in very good agreement with the 2dF- 
GRS £(r) data, especially with the slope at intermediate 
scales (1< r< 10 /i _1 Mpc). Again, we point out that with 
the current data we do not have enough signal at scales 
< 1 /i~ 1 Mpc, so unfortunately we can not assess whether the 
2QZ £(r) steepens at small scales, as seen in the 2dFGRS 
data. 

These results also show that possible deviations from a 
power-law £(r) in the 2QZ are such that they do not become 
evident after being projected along the a direction. 

The determination of f3 assuming the double power-law 
£(r) model can be done as described in Section 4.1. The 
resulting ratio £(s)/£ (r) obtained as a function of separation 
is shown in Fig. |7| By fitting £(s)/£( r ) we find that [3{z — 
1.4) = 0.32± o ;??. The quoted errors are statistical only, and 
do not take into account the systematic uncertainty in the 
(3 measurement due to uncertainty in the form of the £(r) 
model. This is potentially quite large, as demonstrated by 
the difference in the (3 values obtained assuming the single 
and double power-law £(r) models. 

The value of f3(z — 1.4) found is in agreement with 
the values computed from previous estimates, with different 
methods and also based on QSO samples taken from the 
2QZ survey (Hoyle et al. 2002; Outram et al. 2004). The 
value for the linear bias is b(z = 1.4) = 2.84 
a flat fi^j =0.3 cosmology. 



1 49 

; 57 , assuming 



Taking the real- and z-space clustering properties, the 
effect of f3 is to change the amplitude of the correlation func- 
tion at large scales. If, on top of this effect, we add the distor- 
tions due to the small-scales peculiar motions and z-errors, 
we find that the predicted £(s) from our model is a good 
description of the measured £(s) (Fig.|HJ. The reduced Xmin 
is 0.95 (8 d.o.f.), thus indicating a very reasonable fit. 
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Figure 6. The top plot shows the best fitting two power-law 
£(r) model to the w p (a) data, assuming the A cosmology (solid 
line). For scales r < 10 h~ 1 Mpc, the power-law has the form: 
£(r) = (r/6.0)~ 145 , while for scales r > 10ri _1 Mpc, it has the 
form £(r) = (r/7.25) — 2 3 °. The dashed lines show the extrapo- 
lation of the two power-laws to small and large scales. On the 
bottom, the plot shows w p (o)/o as measured from the data (cir- 
cles), overplotted on the projection of the £(r) model from the 
top panel (solid line). 
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Figure 7. £(s)/£(r) measured from the 2QZ survey, taking the 
measured values and the two power-law £(r) model we de- 
rived. As previously, the fit was made on scales where the z-spacc 
distortions are only affected by the large-scale infall and not by 
random peculiar motions of the QSOs. The dashed line and the 
shaded region represent this best fit and its 1 a confidence level. 
This fit corresponds to a value of /3(z = 1.4) = 0.32 
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Figure 8. The plot shows the £(s) results from the 2QZ overplot- 
ted with the £(s) predicted by the double power-law £(r) model, 
after adding the z-space distortions due to peculiar motions. To 
quantify these, we took values for < w 2 > x / 2 and (3(z = 1.4) of 
SOOkms -1 and 0.32, respectively. 



5 COMPARISON WITH RESULTS FROM 
OTHER SURVEYS 

The high statistical significance of the 2dFGRS sample 
makes it unique, and an excellent samp le with which to com- 
pare our results. lHawkins et al.l J2003I) measured the z-space 
correlation function £(s) and the projected correlation func- 
tion, w p (a). Unlike the case of the 2QZ, direct inversion of 
w p (a) of the 2dFGRS allowed measuring the real-space cor- 



relation function f(r) with great accuracy, up to scales of 
~ 20/i _1 Mpc. f(s) is very well defined by a double power- 
law model, up to scales ~ 20 h~ Mpc, where the data stops 
following the two power-law fit and rapidly tends towards 
zero. The flatter behaviour at small scales (< 4 h~ 1 Mpc) 
can be explained by the effects of peculiar velocities with 
dispersion ~ 500 kms -1 . From £(s)/f(r) and £(<r,7r) fitting, 
a value of f}(z) = 0.491^09 was found. 

Although the double power-law model was motivated 
by the form of the 2dFGRS £(»"), we have proceeded by fit- 
ting the model to the 2QZ w p (a) and £(s). Therefore, it is 
now worthwhile checking how the double-power law model 
fitted the 2QZ data compares to the 2dFGRS correlation 
functions. Of course, the form of the high-z QSO correlation 
function need not be consistent with the low-z galaxy corre- 
lation function. However, assuming that the bias is scale in- 
dependent and that the range of scales fitted are not affected 
by non-linear effects, it should be expected the form of £(r) 
to be the same. We compare the measured 2dFGRS w p (cr) 
and £(s) with the predictions from the f (r) 2QZ model. This 
is represented in Fig.|U] 

The shaded region, in both plots, represents the la er- 
rors from the results presented in lHawkins et all i2003l) . The 
diamonds are the measured w p (a)/a and from the 2QZ. 
The dashed and solid lines show the predictions from the sin- 
gle and double power-law £(r) models that best fit the 2QZ 
w p (a) data. The upper plot shows their projection along a 
and the lower plot the respective £(s) functions, with distor- 
tions parametrizes by < w\ > x ' 2 = 800 kms -1 and (3 = 0.87 
(for the dashed line) and (3 = 0.32 (for the solid line). It 
can be seen that the 2dFGRS w p (a) results are very similar 
to the 2QZ and hence to the projection of our £(r) double 
power-law model, even at scales (>20/i _1 Mpc) where the 
2dFGRS d ata does not follow the best-fitting single power- 
law model. ICroom et al.l J2005I) found that the 2dFGRS and 
2QZ samples have the same clustering amplitude. Our re- 
sults for w p (a) corroborate this conclusion. 

There is an offset between the 2QZ £(s) models and the 
2dFGRS data, at small scales. This is due to the different 
z-space distortions in the two data-sets. At large scales, a 
great level of consistency is seen between the results of both 
surveys and the double power-law model. The dotted and 
dash-dotted lines in this plot are the predicted f (s) functions 
for the 2dFGRS, taking the two £ (r) models derived for the 
2QZ (the single and double power-law f(r) models respec- 
tively) and adding the distortions qu antified by the values o f 
0(z = 0.11) and < w\ > 1/2 found bv lHawkins et all fcOOSft . 

The dotted line is not a good representation of the 2dF- 
GRS data since, at 6<s<20ft _1 Mpc scales, as it underes- 
timates the 2dFGRS clustering, considering the 1<t errors. 
The dash-dotted line, that corresponds to the 2dFGRS £(s) 
prediction assuming the double power-law model, is a bet- 
ter description of the clustering of the 2dF galaxies at large 
scales, although it underestimates the clustering at scales 
<3 ft" 1 Mpc. 



COMPARISON WITH CDM MODEL 
PREDICTIONS 



iPeacock fc Doddsl <|l996fl presented a model to describe the 
evolution of the power-spectra (P(k)) of density fluctuations 
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Figure 9. The top plot shows the w p (a) results for the 2dFGRS 
and 2QZ surveys, together with the predictions from the single 
and double power-law models previously derived for the 2QZ. 
The shaded region represents the la error margin for the 2dFGRS 
w p (cr)/cr. The diamonds are the w p (a)/a values measured for the 
2QZ. Similarly, the bottom plot shows the £(s) results, where the 
shaded region is again the la error margin for the 2dFGRS results 
and the diamonds the 2QZ data. Refer to the text for a complete 
description of the models here represented. 



in the non-l i near r egime. This work was then followed by 
ISmith et alJ <l2003f) . who derived a more accurate method 
of including non-linear effects in the power-spectrum, which 
can be applied to more general power spectra. This new 
method of modelling the power-spectrum is based on a fu- 
sion of the halo model with the scaling between lin ear and 
non- linear scales proposed bv lHamilton et all il99ll) . Smith 
and collaborators also presented a set of high-resolution N- 
body simulations that they used to test their mass power- 
spectrum formula. Their results suggest that the model for 
the mass power-spectrum is a good description of its clus- 
tering in a very wide range of scales, k ~ 0.1 — 100 /iMpc - . 

Here we compute the real-space non-linear power- 
spectrum for the mass using the method described in 
ISmith et alJ (I2003T) . By scaling the obtained power-spectrum 
with a linear bias b, we can obtain the predicted CDM 
power-spectrum for the 2QZ catalogue. We also compute 
the z-space power-spectrum by a dding the z-space distor- 
tions to the model, as described in lPadilla fc Bauehl i2002t) . 
After computing the real- and z-space power-spectrum, the 
respective correlation functions can be derived by a simple 
Fourier transform of the P(k) output. 

The top plot in Fig. 1101 shows the £(r) prediction de- 
rived from a ACDM power-sp ectrum, computed using the 
method of ISmith et al.l <l200&t) (dashed line) with parame- 
ters: n° m = 0.3, Q.% = 0.7, as = 0.85, F = 0.17, b = 2.3, < 
w\ y 1 / 2 — 800 kmT 1 and z = 1.4. The solid line is our dou- 
ble power-law model that we derive solely based on the 2QZ 
£(s) and w p (a) measurements. The dotted line is the single 
power-law model derived for the 2QZ. As it can be seen, 
there are reasonable similarities between the ACDM predic- 
tion and the double power-law model. The slopes of the two 
power-laws that describe £(r), fitted from the w p (a) data, 
are close to the predictions from ACDM. The single power- 
law model does not appear to be as close to the ACDM pre- 
diction. The bottom plot shows the £(s) results, after the 
distortions due to j3 — 0.32 and < w\ > 1 / 2 = 800 kms -1 are 
included. Again, the dashed line represents the ACDM P(k) 
prediction. The solid line is the £(s) derived from our double 
power-law f(r) model, and the dotted line the £(s) derived 
from the single power-law model, once the z-space distor- 
tions are added. The circles represent the f (s) measured di- 
rectly from the 2QZ, and the corresponding error bars. The 
value of Xrnin between the £(s) ACDM prediction and the 
£(s) data is 0.87 (8 d.o.f.). Thus, we can conclude that the 
ACDM clustering predictions provide a good description of 
the data, and that the ACDM real-space clustering shows 
more resemblance to the double power-law model than to 
the single power-law model for the QSO £(V). 



7 PARAMETER CONSTRAINTS FROM 
Z-SPACE DISTORTIONS IN £{p,it) 

In Section 4 we estimated /3(z = 1.4) from a direct compari- 
son between f (s) and the £(r) models derived from the w p (a) 
data. We concluded that the /3(z = 1.4) value is strongly 
dependent on the £(r) model assumed. The modelling of z- 
space distortions and their use for parameter constraints can 
be taken a step further, based on the 2-D z-s pace correla- 
tion fu nction £(cr, tt), developing the method of iHovle et alJ 
( 2002) used for the 2QZ lOfc catalogue. 
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Figure 10. The top plot shows the derived biased §(r) model 
from a ACDM P(k) prediction (dashed line), overplotted on the 
best fitting double and single power-law £(r) model fits to our 
w p ((j) data (solid line and dashed line, respectively). There is 
good agreement between the solid and dashed lines. In the bottom 
plot are the respective £(s) functions from the models above, when 
the distortions parametrizes by < wl >V 2 = 800kms~ 1 and = 
0.32 (for the double power-law model) or /3 = 0.87 (for the single 
power-law model) are added. The dashed, solid and dotted lines 
represent the same models as in upper plot, the difference being 
that now the z-space distortions are included. The circles are the 
measured f (s) from the 2QZ data. The ACDM P(k) predicted 
£(s) shows good agreement with the data, considering the errors. 
Approximations made when including the distortions in the two 
power-law model cause a small offset between the solid and dashed 
lines. 



QSO peculiar velocities lead to distortions in the £ (a, n) 
shape. At small scales in a, the random peculiar motions of 
the QSOs cause an elongation of the clustering signal along 
the 7T direction. The predominant effect at large scales is the 
coherent infall that causes a flattening of the £(cr, n) contours 
along the tt direction and some elongation along a. 

Geometric distortions also occur if the cosmology as- 
sumed to convert the observed QSO redshifts into distances 
is not the same as the true, underlying cosmology of the 
Universe. The reason is because the cosmology dependence 
of the separations along the redshift direction is not the 
same as f or the separations measure d in the perpendicular 
direction jAlcock fc Paczvnsdll97gl) . 

In the 2QZ sample it was found that, at small scales, 
there is not enough signal in the data to constrain the ve- 
locity dispersion of the 2dF QSOs with this method. Nev- 
ertheless, it should be possible to determine constraints not 
only on @{z), that parametrizes the large-scale infall of the 
QSOs, but also on the value of Q.^. 

7.1 The 2QZ 7r) 

The measurement of £(cr, tt) is done the same way as £(s), ex- 
cept that now the number of pairs is binned in two variables, 
rather than one. In Fig. llll we display £(<r, tt) determined for 
the 2QZ catalogue. 

The £(<t, tt) contours appear elongated along tt, which 
can be attributed to the effects of peculiar velocities. The 
data is also quite noisy, particularly at small scales in a. The 
errors are computed using the Poisson estimate. It is there- 
fore easy to understand the reason for the noisier contours, 
at those scales: since the a direction actually represents two 
dimensions, while tt only one, the volume of each bin, and 
hence the number of pairs increases as the separation a in- 
creases, leading to a decrease in the Poisson error. 

7.2 Comparison Between Two z-Space Distortions 
Models 

To understand and quantify the effects that are actually 
shaping the £ (a, tt) contours, we develop models for the z- 
space distortions and compare them to the observations. 
There are different models to explain and describe z-space 
distortions in £(cr, tt), based on different assumptions and ap- 
proximations. However, one key assumption in all of them 
is the non-scale-dependence of the bias. 

Here, two models are described and compared. The ba- 
sis for the construction of the two models is the same: we 
start with a model for the real-space correlation function 
£(r), include the effects of the coherent large-scale infall, 
and then convolve the result with the distribution of the 
small-scale pairwise peculiar velocities. 

7.3 z-Space Distortions - Model I 

The large-scale coherent infall of the QSOs, is described in 
Fourier space by (Kaiser 1987; Hawkins et al. 2003): 

P,(k) = (l+0(z)nl)P r (k), (11) 

where P s (k) and P r (k) are the power-spectrum in redshift 
and real-space, respectively, and fik is the cosine of the angle 
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Figure 11. f(<r,7r) measured with the Hamilton estimator. On 
top is the result obtained if an EdS cosmology is assumed, on the 
bottom if the A one is assumed, instead. 



between the wavevector k and the line-of-sight. Translated 
to real-space, these results take the form (Hamilton 1992; 
Matsubara & Suto 1996): 



+ |/J( Z ) + i/3(-) 2 ) Co (r)Po(AO (12) 
(|/3( 2 ) + ^) 2 )6(r)P 2 (M) (13) 



35 



(14) 



where p is now the cosine of the angle between r and it and 
P;(/i) are the Legendre polynomials of order I. £o(r), &(r) 
and £4(7") are the monopole, quadrupole and hexadecapole 
components of the linear £(r) and their form will depend 
on the £(r) model adopt ed. In general, they are given by 
jMatsubara fc Sutdll996T) : 



xdx 



(15) 



The £(<r, 7r) model is then convolved with the pairwise 
peculiar velocity distribution to include the small scale z- 
space effects due to the random motions of the QSOs. Here 
we assume t hat this can b e well described by a Gaussian 
distribution ( Ratcliffelll996T): 



1 



■ exp 



2tt < w'i >V2 " ' \ 2 < wl > 
Hence, £(<r, n) is then given by: 



£(0-,ff)= / £,'{cT,K -W z (l + Z )/H{z))f{w z )dw z 



(16) 



(17) 



where f'(cr, n — w z (l + z)/H(z)) and f(w z ) are given by 
equations 12-14 and 1161 



7.4 z-Space Distortions: Model II 

£(cr,7r) can be defined as (Peebles 1980 ; Hoyle 2000): 

poo 

1 + Z(a,v) = {l+t{r))f{w z )dw z (18) 



Here, it is assumed that the pairwise peculiar velocity 
distribution f(w z ) is a slowly varying function with r. The 
form of f(w z ) is the same as in Model I and is given on 
equation ED The effects of the bulk motions can be included 
in the following way: 



l + e(<r,7r) 



(1 + £(r))/(ui,(l +*)- v(r z ))dw z , (19) 



where v(r z ) is the model used for the bulk motions, as a 
function of t he real-space separa tion along the 7t di rection - 
r z . Following lHale- Suttonl il990T) and lHovld feOOOD . 



v{r z ) = -- 



-n m (z) 0(i H(z)r z 



(20) 



£(cr, 7r) computed with these two models is shown in 
Fig. 1121 The solid line refers to Model I and the dashed 
line to Model II. The input £(r) form is the same for both 
models (£(r) = (r/5.0) -1 ' 8 ). For the plot on the right the 
model is computed considering < u>1 > 1 ' 2 = Okms -1 and 
(3{z) = 0.4. The overall effect is the same for the two 
models, and consists of a compression of the £(cr, 7r) con- 
tours along the n direction and a small elongation along 
the a direction. The plot on the left displays the case with 
< w 2 z > 1/2 = 800 km s" 1 and f3(z) = 0.4. One can see that, 
at small scales, the distortions caused by the velocity dis- 
persion dominate, while at large scales, the distortions due 
to a non-zero (3(z) dominate. The adopted redshift for these 
models is z — 1.4, which is the median redshift of the 2QZ 



© 2004 RAS, MNRAS OOO.IiTITTI 



12 J. da Angela et al. 




Figure 12. £(<j, w) computed with Model I (solid line) and Model II (dashed line). The plot on the left shows the difference between 
the two models when the distortions are quantified by < w 2 > 1 / 2 = Okms -1 and 8(z) = 0.4. The plot on the right also accounts for the 
distortions caused by a non-zero velocity dispersion, and the values < w\ >V 2 = 800 kms -1 and P(z) = 0.4 were used. The £(r) model 
is in all cases a simple power-law of the form (r/5.0) -1 ' 8 . 



survey. These plots also show significant differences in the 
shape of the distortions caused by j3{z) between the two 
models. 



7.5 Are Models I and II self-consistent? 

One simple test can be performed to check if a given model 
is self-consistent. Averaging £(cr, ir) in several annuli will 
give £(s). In addition, at large scales, z-space distortions 
are mainly affected by the large-scale coherent infall, which 
implies the relation between £(s) and £(r) given by equa- 
tion EI] Hence, if a model is self-consistent, after averaging 
£(<7, 7t) at constant s, the final result should be the same as 
the initial input £(r), scaled with (1 + 2/3(3(z) + l/5/3(z) 2 ). 
If not, then the model must be wrong. Here, this test is done 
for Models I and II. 

Fig. 1131 shows the result of this test. The circles re- 
fer to Model I and the diamonds to Model II. In both 
cases a ^(a, n) model is computed using an input power-law 
£(r) = (r/5.0)" 1 - 8 , P(z) = 0.4 and < w\ > 1/2 = 800 kms" 1 . 
Then, the average £(cr, it) is computed in constant annuli 
of s, thus obtaining £(s). The output £(s) is then divided 
by the predicted £(s) using equation 1101 One can see that 
Model I seems self-consistent, as the predicted £(s) from £(r) 
matches with the linear scaling of £(r), through {3(z), at large 
scales. However, this is not the case for Model II. There is a 
discrepancy between the predicted £(s) and the one derived 
from averaging the z-space distortions on £(er, n). This offset 
corresponds to an offset on the value of (3(z) of 0.25, which 
is quite significant. Given this result, only Model I will be 
considered, hereafter. 
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Figure 13. §(s) computed from averaging §(<r, it) in annuli of 
constant s = \/o" 2 + 7r 2 and divided by the power-law predicted at 
large-scales. The circles correspond to Model I and the diamonds 
to Model II. It was assumed f(r) = (r/5.0) -1 - 8 , /3(z) = 0.4 and 
< ui 2 >!/ 2 = 800 kms -1 . £(s) computed from Model I seems to 
agree very well with the prediction from the linear regime scaling 
from £(r), given by equation 1101 The same does not happen with 
Model II, where an offset between the predicted and the derived 
£(s) is observed. It can also be seen for both models the effect of 
the velocity dispersion, that causes a flattening of £(s) at small 
scales. 



7.6 Including Geometric Distortions 

It is useful to make some definitions before describing the 
cosmology f i tting through geometric distortions. Following 
iHovle et alJ (I2002T) . let the true, underlying cosmology of the 
Universe be the true cosmology, the cosmology used to build 
the model £(cr, 7r) the test cosmology, and the cosmology as- 
sumed to derive the r — z relation, both in the model and 
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Figure 14. £(cr, 7r) computed using Model I for different assumed and test cosmologies. The left plot shows the case where an EdS 
cosmology was assumed, to measure the correlation function, but the underlying cosmology of the Universe has Q,^ = 0.3 and Sl^ = 0.7. 
On the plot on the right is the case where the assumed cosmology is a flat cosmology with Q ^ = 0.3 but the underlying cosmology 
has £7^ = 1-0- Assuming values for Qj^ higher /smaller then the true ones, will cause a flattening/elongation of the £(cr, n) along the 
line-of-sight, besides the overall effect on the amplitude of the correlation function, also visible in these plots. In these cases, z-spacc 
distortions due to peculiar velocities were disregarded. 



the data to measure the correlation function, the assumed 
cosmology. 

Since we are comparing the geometric distortions in 
both the data and the model relative to the same assumed 
cosmology, the test cosmology used in the model that best 
matches the data should be the true cosmology of the Uni- 
verse, in the absence of noise. 

In this paper we assume a spatially fiat cosmology, and 
choose to fit the variable Q^i, hence fixing Q% = 1 — Q^. 

The relation between the separations a and n 
in the test and assumed cosmologies (referred to by 
t he subscripts t and a, respectiv el y) is the following 



y) is tne K 
llHovlell200d) : 



( Ballin ger. Peacock fc Heavenslll99a) . 
crt — f±a a = —o~a (21) 

Ha 

At 

7T( = /||7T a = — 7T a (22) 

where A and B are defined as follows (for spatially flat cos- 
mologies) : 



.4 



Ho ^ + ^(1 + ^)3 

dz 



HoJ + fi° I (l + z') 3 



(23) 



(24) 



In the linear regime, the correlation function in the as- 
sumed cosmology will be the same as the correlation function 
in the test cosmology, given that the separations are scaled 
appropriately, i.e.: 



it{o~t, 7T t ) = fa(Cai 7Ta). 



(25) 



In Fig. I14l the geometric distortions due to assuming a 
cosmology different from the true, underlying cosmology are 
represented. On the left plot is £(cr, 71") derived from Model I 
using an f2° n =0.3 test cosmology, while an EdS cosmology 
was assumed to derive the comoving distances. The plot on 
the right shows the case where the underlying cosmology has 
fl^ = 1.0 but a A cosmology was assumed. So, assuming a 
value of fi^, higher than the true one will cause a compres- 
sion of the £(cr, 7r) contours along the line of sight, whilst if a 
too small value for £7°^ is assumed, it will cause an elongation 
of the contours along the line of sight and a high clustering 
amplitude to be observed. Redshift-space distortions due to 
the effects of peculiar velocities were not included in these 
models. 

Ignoring for the moment the distortions caused by the 
peculiar velocities and the effects of noise, the observed 
shapes of the measured and modelled £(<r, n) will be the 
same when the test cosmology matches the true, underly- 
ing cosmology of the Universe (for whatever cosmology is 
assumed to convert the redshifts to comoving distances). 
Hence, by fitting the geometric distortions in the data, one 
should be able to determine the true cosmology. 

7.7 Fitting Procedure 

The fitting proced ure was developed from the one used by 
iHovle et alJ <|2002T) . In summary, for a given value of /3(z), a 
£(<7, 7r) model is generated in a chosen test cosmology. Then, 
the separations a and n are scaled to the same cosmology 
that was assumed to measure the actual data. The final 
model for £(cr, 7r) is then compared to the data. This method 
is repeated for different test cosmologies and values of /3(z). 
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Since we are fitting distortions in the shape of the z- 
space correlation function, the correct spherically-averaged 
amplitude of £(r) must be given as an input to the model. 
Otherwise the fit will be driven by offsets in the amplitude 
of the £(<t, 7r) model from the data, rather than the shape 
distortions, which would introduce systematic errors in the 
constraints obtained for and (3{z). 

The following steps are taken in the fitting procedure: 

• Assume a cosmology and measure £(s), w p (a), £(a,Tr). 

• Take a model for the real-space correlation function, e.g. 
a double power-law. This model should be a good description 
of the observed data, £(s) and w p (a). 

• Choose a pair of test values of and f3{z). 

• The model for £(r) is a good description for the data 
in the assumed cosmology. What is actually needed at this 
stage is a £ (a, n) model in some test cosmology, hence the 
correct input for this model is £(r) in that same test cosmol- 
ogy. Since, in the linear regime, ft(r t ) = £ a (r a ), one has only 
to compute the real space separation in the assumed cosmol- 
ogy to get £t (r t ). r t is given by r t = y 'of + (n t - w z /H t ) 2 
and the relation between r a and r t is: r a = r t /(fj_f\\) 1 '' 3 . 

• Using that model for £t(r t ), compute ^ t {a t ,-K t ). Then, 
include the geometric distortions by scaling (,t{o't,'^t) back 
to the assumed cosmology, in a similar way as described in 
the previous step. To get £ a (o"a, n a ), one needs to scale the 
separations a t and n t to a a and n a , using equations 1211 and 

m 

• Adding the effects of large-scale infall not only intro- 
duces distortions in £(cr, it) but also shifts the amplitude of 
the correlation function, by an amount that depends on the 
value of /3(z) taken. Since the amplitude of the spherical- 
averaged correlation function must remain the same (i.e. 
match the £(s) data), whatever /3 and are used as test 
values, the amplitude of the input £(r) model is allowed to 
vary in the fit, guaranteeing that the fit is being made to the 
distortions in £(cr, n), since the averaged amplitude remains 
the same for whatever combination of /3 and 

• For the best fitting value of this amplitude factor, deter- 
mine the x 2 value for the fit of this model to the data. In this 
step we are assuming that the £(<r, it) bins are independent. 

• Repeat this procedure for different combinations of 
and f3(z). 

The number of degrees of freedom in the \ 2 fit is the 
total number of bins where £(<r, tt) from the model is fitted 
to the data minus the number of free parameters. If the fit 
is to fi^j and f3(z) and the averaged amplitude of 7r) 
is allowed to float so it matches £(s), the number of free 
parameters will be three. The velocity dispersion was fixed in 
these fits: taking into account the z-errors of the survey and 
the intrinsic velocity dis persion of the QSOs , we assumed 
< wl > 1 ' 2 = 800 kms" 1 llOutram et alJl2003D . 



7.8 Further Constraints on and /3(z) from 
QSO Clustering Evolution 

All the other factors being the same, the greater the value 
of the true, underlying f2„, the more elongated the £(tr, 7r) 
will be along 7r; and the greater the value of [3(z), the flatter 
these contours will be. Hence, a degeneracy is expected to 
occur in the confidence levels in the (3{z)\ plane. In 



order to get better constraints on fi^n an d also f3(z), this 
degeneracy needs to be broken. 

A possible way to break the degeneracy is to combine 
these results with a constraint derived from consideration 
of QSO clustering evolution. From the value of the mass 
correlation function, at z = 0, linear perturbation theory 
can be used, in a given test cosmology, to compute its value 
at z = 1.4. Then, considering that the bias (b(z = 1.4)) is 
given by the ratio of the QSO £(r) and the mass £(r), and 
that the former is, at large scales, related to the measured 
QSO £(s) through f3(z — 1.4), estimates of this parameter 
can be obtained for a given test cosmology. 

To compute the mass correlation function amplitude at 
low z, we can use the values of £(s) and found for the 
2dFG RS survey. Using the values found bv lHawkins et all 
(2003), the first step is to compute the respective value of 
b(z = 0). 

For each test cosmology, the bias parameter of the 
galaxies at z = is: 



b(z = 0) = 



(26) 



0{z = 0) 

Consider now the volume averaged two-point correla- 
tion function £ s given by: 



f* 4ns' 2 ds' 



(27) 



For the 2dFGRS, £(s ) is found to be well d escribed by a 
double power-law model (Hawkins et al. 2003). To compute 
equation 1271 for the 2dFGRS, that model can be used in the 
integral on the numerator. Non-linear effects due to peculiar 
velocities in the sample should be insignificant by taking the 
upper limit of the integral s — 20 h~ Mpc. 

Then, the equivalent averaged correlation function in 
real-space can be determined by: 



?{z = 0) = 



1 + U{z = 0) + U{z = 0) 2 



(28) 



Now the real-space mass correlation function is ob- 
tained with: 



£,mass ( z 0) 



r> = o) 

b(z = 0) 2 ' 



(29) 



where b(z = 0) is given by equation 1261 

Once determined the real-space correlation function of 
the mass at z = 0, its value at z = 1.4 is obtained using 
linear perturbation theory. Hence, at z = 1.4, the real-space 
correlation function of the mass will be: 



^mass (z 0) 

G(z = 1.4) 2 ' 



(30) 



where G(z) is the growth factor of perturbations, given by 
linear theory an d depends on cosmology (in t his case the 
test cosmology) dCarroll. Press, fc Turnei1ll992l) . 

Once the value of (,lnass{z = 1.4) is obtained for a given 
test cosmology, the process to find fi(z = 1.4) is similar to 
the one used to find ^ a ss( z — 0), but now the steps are 
performed backwards: 

£ s (z — 1.4) can be measured in a similar way as £ s (z — 
0). The bias factor at z m 1.4 is given by: 



b(z = 1.4) 2 = 



r 



1.4) 



1.4) 



(31) 
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where £m ass is given by equation 1301 arid £ r (z = 1.4) is ob- 
tained by: 



Q m assumed = 1 .0 



?{z = 1.4) = 



V{z = lA) 



1+ §£(* = 1.4) + |/3(0 = 1.4)2 



(32) 



The value of (3(z = 1.4) can then be determined by 

(n m (z = i.4)) - 6 



/3(z = 1.4) = 



(33) 



6(2 = 1.4) 

where b(z = 1.4) is given by equation 1311 and Q m (z = 1.4) 
is the value of the matter density at z — 1.4, given by: 



(34) 



n^(i + 2) 3 + n° ' 

for a fiat universe. 

In the end, for a given value of in the test cosmology, 
(3(z) will be obtained by solving a second order polynomial 
equation. 

The confidence levels on the computed values of j3{z = 
1.4) can be obtained by considering the errors on this calcu- 
lation. These are found by identifying the factors that con- 
tribute to the error, and adding the components in quadra- 
ture. Here, the components contributing to the error on 
(5(z = 1.4) are (5(z = 0), £ s (z = 0) and £ s (z = 1.4). 

7.9 Results 

We assume the double power-law £(r) model to compute 
the constraints on fi^ and f3(z) from the f(<r, 7r) shape. By 
substituting this function in Eg. 1151 we find the subsequent 
expressions for the moments of the correlation function and 
hence the form of £(er, it). 

It should be noted that there is some sensitivity in the 
£(u, 7r) fits to the detailed form of the assumed QSO £(r). 
If the model for £(r) is inaccurate then instead of fitting 
the distortion's shape in the [a, n] plane, the fits will be 
dominated by small differences in the average amplitude. As 
has been seen already in Section 4, the (3 constraints from the 
£(s) /£(»") ratio are also sensitive to the assumed £(r) form. In 
both cases, assuming single power-law models for £(r) would 
lead to higher fitted values of f3(z — 1.4). We have argued 
that the single power-law model narrowly fails to fit the 
observed QSO £(s) and w p (a) results. As well as improving 
the fits to these data, the double power- law model for £(r) 
is also a better representation of the correlation function 
results in the 2dFGRS (and Durham/UKST) galaxy survey 
and also in ACDM simulations. 

The solid lines and the shaded areas in Fig. 1151 repre- 
sent the results from fitting £(er, 7r) models to the z-space 
distortions. The plots show similar likelihood contours for 
both the assumed cosmologies. This is as expected, since 
they should in theory be independent of the assumed value 
of The small differences that are seen are due to dif- 
ferences in binning the pairs between the two assumed cos- 
mologies. It can be seen that the constraints on f3(z — 1.4) 
are much stronger than the constraints on f2^. Also, the 
lower the value of O^, i.e. higher the more significant 
the effects of the geometric distortions are. This is reflected 
in the increased curvature of the contours for low f2^. This 
positive slope of these contours helps breaking the degen- 
eracy between and (3, when combining the results with 
constraints from linear growth theory (see below). 
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Figure 15. The confidence levels in the [O^, fi(z)\ plane obtained 
if f(r) is described by the double power-law model, obtained from 
fitting the distortions in the £(<r, 7r) contours (grey-scale and solid 
confidence levels) and from QSO clustering evolution derived from 
linear theory (dashed la and 2a confidence levels). The dotted 
line shows the joint two-parameter la confidence level obtained 
from combining both methods. 



The dashed lines represent the la and 2a levels obtained 
from the bias evolution with redshift. We note that there is 
excellent overlap with the la constraints from z-distortions. 
Essentially, it shows that the amplitudes of mass clustering 
at 2 m 1.4, allowed by 2dFGRS dynamical analysis at z ~ 
0.1, are consistent with the 2-space distortions results from 
the QSOs, also at z w 1.4. 

Finally we combine both the constraints from z- 
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distortions and evolution to obtain the joint constraint rep- 
resented by the dotted line. It shows the joint two-parameter 
la confidence level obtained with both methods. 

If the EdS cosmology is assumed, the best fitting val- 
ues for and /3(z = 1.4) from the joint constraints are: 
Q,° m = 0.20±o:?i and f3(z = 1.4) = 0.70t^8, with a reduced 
Xmin °f 1-10 (12 d.o.f). If the A cosmology is assumed in- 
stead, the best fitting values for and /3(z = 1.4) will 
be: Q,° m = 0.35±{J;" and P( z = L4 ) = °- 50 ±o:i5> with a 
reduced Xmin °f 1-05 (17 d.o.f.). Making the assumption 
of a double power-law model and the A cosmology, the 
P(z) constraint obtained from fitting the shape of £(<x, 7r), 
/3(z = 1.4) = 0.50t° ", is slightly higher than the value of 
P{z = 1.4) = 0.32to.ii obtained from fitting £(s)/f(r), but 
these results are consistent within the margin of error. 

Our results are also consistent with previous estimates 
of (3{z = 1.4) and Q. m from the 2QZ. This can be s een by 
comparing Fig. [T^l with Fig. 9 of IHovle et all J2002ft . How- 
ever, the £( < t, 7r) a nalysis presented here supersedes that of 
IHovle et alj ]2002l) , not only because of the ~ 2 x larger QSO 
data set but also because o f our improved z-spa ce distor- 
tion analysis. More recently lOutram et alJ <l2004l) used the 
power spectrum to study the z-space distortions in the 2QZ 
clustering. Their method probes larger scales than those in 
this work. Their constraints of (3{z = 1.4) = 0.451JJ;™ and 
Oto = 0.29^0 ^9 are also in very good agreement with our 
estimate. 



8 CONCLUSIONS 

To model z-space correlation functions to constrain Q m and 
/3(z = 1.4), we have found that we need an accurate de- 
scription of the real-space QSO correlation function. Al- 
though a simple power-law model gives a reasonable fit 
to w p (a), we find other evidence that the real-space cor- 
relation function, £(r), may show non-power-law features. 
We have argued that fitting w p (a) and £(s) simultaneously 
requires more than just a single power-law moderated by 
the effects of peculiar velocities. Further motivated by the 
non-power-law £(r) shape measured for the 2dFGRS sur- 
vey, we have found that a similar double power-law model 
can explain the observed 2QZ £(s) and w p (a). The parame- 
ters of the best-fitting model are: ro = 

6.0toi /i _1 Mpc and 
7 = lA5t° i 7 r, fOT r < 10/l _1 Mpc, and r = 7.25/t _1 Mpc 
and 7 = 2.30^^3, for r > 10/i _1 Mpc (assuming A). 

We then compared both the single and double power- 
law £(r) models with the 2dFGRS results and theoretical 
predictions from CDM clustering models. When compar- 
ing the 2dFGRS and 2QZ projected correlation functions, 
we found that our proposed double power-law model lies 
within the 2dFGRS errors, which suggests that this 2QZ 
£(r) model provides an acceptable description of the 2dF- 
GRS £(r). With the assumed A cosmology, even the ampli- 
tudes of the 2QZ and 2dFGRS correlation functions are sim- 
ilar; with the assumed EdS cosmology the 2QZ amplitude 
is a factor of two smaller at ~ 10ft _1 Mpc scales. The mass 
correlation function amplitude evolves with red shift, and the 
same is probably true for galaxies and QSOs llCroom et alJ 
120051) . Therefore, the agreement of the amplitudes in the 
A case is probably a coincidence. However, the 2QZ and 



2dFGRS correlation functions are certainly similar in their 
shape. 

The ACDM £(s) prediction appears to be a good fit to 
the 2QZ £(s). We also found that the £(r) double power- law 
model that we derived from the 2QZ clustering results is a 
very good match to ACDM £(r). On the other hand, the 
single power-law 2QZ £(r) model, derived from the w p (a) 
data, produces a poorer match to the 2dFGRS data and the 
ACDM predictions. 

Using the double power-law £(r) model and the £(s) 
measurements, we estimated (3{z — 1.4) = 0.32^q'^, by 
computing and fitting f(s)/£(r). This result contrasts with 
the higher value of f3(z = 1.4) = 0.8Tt{5;jj?, found when we 
consider the single power-law £(r) model. Similar effects are 
seen in 2dFGRS, where the single power-law £(r) model is 
also rejected by the data. 

Modelling the z-space distortions in £(cr, n) allows con- 
straints on f2„ and f3(z = 1.4) to be derived. For the 
analysis we have assumed a pairwise velocity dispersion of 
800 km s~ 1 . Since this is dominated by z-errors (±600 km s~ 1 
pairwise) our results should be robust to reasonable varia- 
tions in this parameter; the pairwise velocity dispersion for 
the mass, for = 0.3 at z = 1.4, is < w't > 1/2 ^ 400 kms -1 
iHovlell2000l) and under simple assumptions this may only 
increase to w 700 kms -1 , for = 1.0. 

We compared the £ (a, n) z-space distor tions model with 
that previously used bv lHovle et al] <l2002l) . and found that 
the latter suffers from a systematic offset between the ampli- 
tude of the input £(s) and the £(ct, 7r) amplitude. This does 
not seem to occur in the model used in the current work 
jMatsubara fc Sutolll99^ . 

Although the £(cr, n) fitting is sensitive to the form of 
the £(r) model, assuming the double power-law model de- 
rived above, we have obtained useful constraints on /3(z = 
1.4) from the z-space distortions analysis. The constraint 
on is not as strong from z-space distortions alone, be- 
ing quite degenerate with the constraints on /3(z = 1.4). 
Combining these z-space distortions results with those from 
QSO bias evolution helps to break this degeneracy and 
provides stronger constraints. The resulting constraints are 
Q,° m = 0.35±° -\l and [3{z = 1.4) = 0.50t°:\l, assuming the A 
cosmology. Assuming the EdS cosmology instead produces 
similar results. Thes e values are in very go od agreement with 
the results found bv lOutram et al.l ]2004l) . 



ACKNOWLEDGMENTS 

The 2dF QSO Redshift Survey (2QZ) was compiled by 
the 2QZ survey team from observations made with the 
2-degree Field on the Anglo-Australian Telescope. JA ac- 
knowledges financial support from FCT/Portugal through 
project POCTI/FNU/43753/2001 and from the European 
Community's Human Potential Program under contract 
HPRN-CT-2002-00316, SISCO. PJO acknowledges the sup- 
port of a PPARC Fellowship. We thank Catarina Lobo, Nel- 
son Padilla and Carlton Baugh for useful comments. 

REFERENCES 

Alcock, C. & Paczynski, B. 1979, Nature, 281, 358 

© 2004 RAS, MNRAS 000.HHT71 



Ballinger, W. E., Peacock, J. A., & Heavens, A. F. 1996, 

MNRAS, 282, 877 
Baugh, C. M. 1996, MNRAS, 280, 267 
Calvao, M. O., de Mello Neto, J. R., & Waga, I. 2002, 

Physical Review Letters, 88, 091302 
Carroll, S. M., Press, W. H., & Turner, E. L. 1992, 

ARA&A, 30, 499 
Croom, S. M., Boyle, B. J., Shanks, T., Smith, R. J., Miller, 

L., Outram, P. J., & Loaring, N. S., Hoyle, F., da Angela, 

J. 2005, MNRAS, 356, 415 
Croom, S. M., Smith, R. J., Boyle, B. J., Shanks, T., Miller, 

L., Outram, P. J., & Loaring, N. S. 2004, MNRAS, 349, 

1397 

Davis, M. & Peebles, P. J. E. 1983, ApJ, 267, 465 
Hale-Sutton, D. 1990, PhD thesis, Univ. Durham 
Hamilton, A. J. S., Kumar, P., Lu, E., & Matthews, A. 

1991, ApJl, 374, LI 
Hamilton, A. J. S. 1992, ApJ, 385, L5 
Hamilton, A. J. S. 1993, ApJ, 417, 19 
Hawkins, E., et al. 2003, MNRAS, 346, 78 
Hoyle, F. 2000, PhD thesis, Univ. Durham 
Hoyle, F., Outram, P. J., Shanks, T., Boyle, B. J., Croom, 

S. M., & Smith, R. J. 2002, MNRAS, 332, 311 
Kaiser, N. 1984, ApJl, 284, L9 
Kaiser, N. 1987, MNRAS, 227, 1 

Kauffmann, C, & Haehnelt, M. 2000, MNRAS, 311, 576 

Lewis, I. J., et al. 2002, MNRAS, 333, 279 

Marsh, T. R., et al. 2002, ASP Conf. Ser. 261: The Physics 

of Cataclysmic Variables and Related Objects, 200 
Matsubara, T. & Suto, Y. 1996, ApJ, 470, LI 
Myers, A. D., Outram, P. J., Shanks, T., Boyle, B. J., 

Croom, S. M., Loaring, N. S., Miller, L., & Smith, R. J. 

2005, MNRAS, 321 
Outram, P. J., Hoyle, F., & Shanks, T. 2001, MNRAS, 321, 

497 

Outram, P. J., Hoyle, F., Shanks, T., Croom, S. M., Boyle, 
B. J., Miller, L., Smith, R. J., & Myers, A. D. 2003, MN- 
RAS, 342, 483 

Outram, P. J., Shanks, T., Boyle, B. J., Croom, S. M., 
Hoyle, F., Loaring, N. S., Miller, L., & Smith, R. J. 2004, 
MNRAS, 348, 745 

Padilla, N. D. & Baugh, C. M. 2002, MNRAS, 329, 431 

Peacock, J. A. & Dodds, S. J. 1996, MNRAS, 280, L19 

Peebles, P. J. E. 1974, ApJl, 189, L51 

Peebles, P. J. E. 1980, The Large Scale Structure of the 

Universe, Princeton University Press, 1980 
Ratcliffe, A. 1996, PhD thesis, Univ. Durham 
Ratcliffe, A., Shanks, T., Parker, Q. A., & Fong, R. 1998, 

MNRAS, 296, 191 
Saunders, W., Rowan- Robinson, M., & Lawrence, A. 1992, 

MNRAS, 258, 134 
Shanks, T. & Boyle, B. J. 1994, MNRAS, 271, 753 
Shanks, T., Hale-Sutton, D., & Boyle, B. J. 1988, IAU 

Symp. 130: Large Scale Structures of the Universe, 130, 

371 

Smith, R. E., et al. 2003, MNRAS, 341, 1311 



z-Space Distortions in the 2QZ Survey 17 



© 2004 RAS, MNRAS 000.1111171 



